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We compute the Coulomb distortion produced by an expanding and highly charged 
fireball on the spectra of low transverse momenta and mid rapidity pions produced in 
central high energy heavy ion collisions. We compare to data on Au+Au at 11. 6 A GeV 
from E866 at the BNL AGS and of Pb+Pb at 158A GeV from NA44 at the CERN 
SPS. We match the fireball expansion velocity with the average transverse momentum of 
protons and find a best fit to the charged pion ratio when the fireball radius is about 10 
fm at freeze-out. This value is common to both AGS and SPS data. 
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An important feature to account for in the analysis of the spectra of secondaries pro- 
duced in the collision of heavy systems is the presence of a large amount of electric charge. 
Due to the long-range nature of the electromagnetic interaction, the spectrum of charged 
particles will be distorted even after freeze-out. For central collisions, this Coulomb effect 
can be more significant when there is strong stopping and the participant charge in the 
central rapidity region is an important fraction of the initial charge. 

Another feature to consider is that the field-producing charge distribution is in general 
not static but rather participates in the dynamics responsible for matter expansion after 
the collision. The combined role played by Coulomb distortions and expansion in the 
description of charged particle spectra has been the subject of some recent work. Koch |I[ 
and Barz, Bondorf, Gaardh^je and Heiselberg have developed approximate models to 
describe the situation in which expansion takes place predominantly along the collision 
axis. In this work, we focus on the description of Coulomb effects on pion spectra from 
a spherically symmetric expanding source. We compare our calculation to mid-rapidity 
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pions produced in central Au+Au reactions at 11. 6A GeV from E866 at the BNL AGS 
and in central Pb+Pb reactions at 158A GeV from NA44 at the CERN SPS. A detailed 
analysis can be found in Refs. [[|. 

Before proceeding on to the model calculation, let us say a few words about the as- 
sumption of spherical symmetry. It has been known for some time that the momentum 
distributions of secondaries are somewhat forward-backward peaked, specially at the SPS, 
even for central collisions |4[], and this observation can cast some doubt about the validity 
of a model that assumes a spherically symmetric fireball. However, spherical geometry is 
not essential to the basic physics and can be relaxed at the expense of additional comput- 
ing time. Nevertheless, let us provide the following arguments in favor of its use. First, 
we will be comparing with the transverse momentum distributions at mid-rapidity where 
the impact of spherical asymmetry should be less important than near the fragmentation 
regions. Second, pion interferometry of central Au+Au collisions at the AGS J5] and of 
central Pb+Pb collisions at the SPS both yield comparable values for the transverse 
and longitudinal radii at the time of pion freeze-out. These radii are about twice the 
radius of a cold gold or lead nucleus. Third, as we will later see, the transverse surface 
of the fireball needs to expand with a speed more than 90% that of light in order to 
reproduce the average proton transverse momentum. Since the longitudinal surface of 
the fireball cannot travel faster than the speed of light, this means that in velocity space 
the fireball is nearly symmetric. These phenomena, although not yet measured at that 
time, were already known to Landau 0]. The essential insight from his model is not the 
degree of stopping but rather the point that significant transverse expansion sets in after 
the longitudinal and transverse radii become comparable in magnitude. Thereafter, from 
the point of view of a distant observer, the expansion is not as asymmetric as one might 
originally think. This model was later developed by others, including Cooper, Frye and 
Schonberg ||. 

A uniformly charged sphere which has a total charge Ze and whose radius R increases 
linearly with time t from a value R at time t at a constant surface speed v s produces 
an electric potential 



In the center-of-mass frame of the fireball the charge moves radially outwards, hence 
there is no preferred direction and consequently the magnetic field produced by this 
moving charge configuration vanishes. The fireball parameters are related by Rq = v s to. 
If / ± (r, p,t) represents the ±e test particle phase space distribution then, when ignoring 
particle collisions after decoupling, its dynamics is governed by Vlasov's equation. 



where E p = ^p' 2 + m 2 , m is the meson's mass and E(r, t) = — \7 r V(r, t) is the time- 
dependent electric field corresponding to the potential V(r, t). 

The solution is found by the method of characteristics. This involves solving the classical 
equations of motion and using the solutions to evolve the initial distribution, taken to be 




r > R = v s t 
r < R = vd 



(1) 




|-.V r ±eE(r,t).V p / ± (r,p,t) = 0, 



(2) 
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thermal 

/ ± (r,p,* ) = e X p{-(£ p ±V(r,* ))/T} , (3) 

forward in time. The pion's asymptotic momentum is computed numerically by a sixth 
order Runge-Kutta method with adaptive step sizes from a set of initial phase-space 
positions. The final momentum distribution is the result of computing the trajectories 
for many initial phase space points. The initial radial position was incremented in N r =50 
steps with spacing Ar/i? = 0.02. The initial momentum was incremented in N p =600 
steps with spacing Ap = 1 MeV. The cosine of the angle between the initial position and 
momentum vectors was incremented in N 2 =100 steps of size Acos# = 0.02. Hence the 
total number of trajectories computed was 3 x 10 6 for each set of initial conditions. 

To determine the value of the surface expansion velocity, we match the average trans- 
verse momentum, as inferred from our assumption of a uniformly expanding sphere, with 
the measured one. The result is 



2 -(2 + ^)0 



v 2 



(Pt)=tt — m P . (4) 

Note that v s should notbe interpreted as a hydrodynamic flow velocity, rather, it embodies 
the combined effects of hydrodynamic flow and thermal motion of the net charge carriers, 
mainly protons. 

For central Pb+Pb collisions at 158 A GeV ||, the number of protons per unit rapidity 
in the central region is on the order of 30. If we consider that the rapidity region spanned 
by the fireball is between 1 and 5 units, then we take as the effective fireball's charge 
Z = 120. To model the primordial distribution, we use an exponential parametriza- 
tion of dN/p 2 dp from which we extract an effective temperature T e // = 110 MeV up to 
rrtrp — m = 500 MeV. Figure 1 shows this representation in comparison to the transverse 
mass distributions of positive and negative pions. Both NA44 and NA49 ||T0| , [lT| have 
reported the transverse mass distribution in central Pb+Pb collisions at mid-rapidity to 
be dN/p T dp T oc exp(— m T /T P ) with T P = 290 MeV. This corresponds to an average 
transverse momentum of 825 MeV. The value v s = 0.916 gives, according to Eq. (|4]), the 
match of the model to the proton spectra. The best fit to the n + /n~ ratio is obtained for 
a value of the freeze-out radius for pions of 10 fm with a Xp.d.o.f. = 1-67. 



For central Au+Au collisions at 11.6^4 GeV [12], the number of participating protons 
may be estimated from the data to be Z = 116. The average proton transverse mo- 
mentum is 820 MeV yielding the surface speed v s = 0.914. The slopes of the 7r + and 
7r~ distributions up to — m = 400 MeV are the same to within several MeV; their 
average is 111 MeV. All of these quantities are remarkably similar to those in central 
Pb+Pb collisions at the much higher energy of the SPS. The best fit to the n~ /ir + ratio is 
obtained for a value of the freeze-out radius for pions of 10 fm with a Xp.d.o.f. = 0.53. The 
computed ratio is compared to the data in Fig. 2 and its appearance is quite satisfactory. 

In conclusion, we have shown that the suppression of the ratio tt + /tt~ in central Pb+Pb 
collisions at the SPS, as observed by NA44, and in central Au+Au collisions at the AGS, 
as observed by E866, can be quantitatively understood as a Coulomb effect generated by 
the electric field of an expanding and highly charged fireball. This ratio provides a good 
measure of the size of the fireball at decoupling. In principle, a different parametrization 
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Figure 1. The transverse mass distri- 
butions at mid-rapidity for 7r + and 7r~. 
The solid lines are exponentials with an 
inverse slope of 110 MeV. The data are 
from NA44 |K|. 
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Figure 2. The ratio 7t~/tt + vs. m^ — mat 
the AGS. The data is from Ref. JT|. The 
fireball radius at pion freeze-out is 10 fm. 



of the primordial distribution, such as a two temperature fit, might lead to an even better 
representation of data. 
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